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Abstract 

The near-threshold np — > cross section is calculated in chiral perturbation theory to next- 
to-leading order in the expansion parameter ^Mm 7r /A x . At this order irreducible pion loops 
contribute to the relevant pion-production operator. While their contribution to this operator 
is finite, considering initial- and final-state distortions produces a linear divergence in its matrix 
elements. We renormalize this divergence by introducing a counterterm, whose value we choose 
in order to reproduce the threshold np — > d7T° cross section measured at TRIUMF. The energy- 
dependence of this cross section is then predicted in chiral perturbation theory, being determined by 
the production of p-wave pions, and also by energy dependence in the amplitude for the production 
of s-wave pions. With an appropriate choice of the counterterm, the chiral prediction for this energy 
dependence converges well. 
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I. INTRODUCTION 



Ever since the days of current algebra, soft-pion theorems, which relate processes in- 
volving different numbers of pions, have been used to understand threshold pion-production 
reactions . These relations originate from the chiral symmetry of the underlying La- 
grangian of quantum chromodynamics (QCD) and have now been systematized into an 
effective field theory called chiral perturbation theory (%PT) jsl Q]. In %PT, there is not 
only an understanding of these relations between amplitudes, but also an organizing prin- 
ciple (power counting), which suggests a hierarchy of contributions to the amplitudes and 
thus makes it possible to assign a theoretical error to a calculation. ;\PT can be regarded 
as an extension of and improvement on current algebra and it has found many applications. 
Significant progress has been made in the pion-pion and pion-nucleon sectors (for a review 
of the latter see and for nucleon-nucleon scattering P, 0, Qj. Here we will apply it to 
threshold pion production in few-nucleon systems. 

The application of chiral perturbation theory to threshold pion-production reactions has 
met with mixed success. In the charged-pion production reactions the leading pion rescat- 
tering mechanism involves the Weinberg- Tomozawa (WT) ttN amplitude. Considering this, 
plus certain higher-order effects, produces moderate agreement between theory and data 
cross sections |8|. Approaches to these reactions based on meson-exchange models are also 
quite successful jj], [l0( . 

Meanwhile, both meson-exchange models and %PT initially struggled to describe the 
pp — > ppir cross section measured close to threshold. Data from the Indiana University 
Cyclotron Facility 11] and the The Svedberg Laboratory in Uppsala, Sweden were five 
times larger than phenomenological calculations [l^ . Two early proposals to remedy this 
deficiency were short-range interactions Q, Q] and irN amplitudes with a particular off- 
shell dependence Q|. In %PT these effects are both higher order, and they are related to 
each other through field redefinitions jl^] • This sensitivity to higher-order effects is partly a 
consequence of the particular isospin structure of the pp — > ppn reaction, which forbids the 
pion rescattering via the WT term that is the leading effect in the charged-pion channels. 
(In %PT there is also a partial cancellation between higher-order one-body and rescattering 



terms, which is not present in meson-exchange model treatments 
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19,t20l.) But 



even at a formal level, it quickly became apparent that the conventional xPT expansion of 
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amplitudes in jj- (with M and m, the nucleon and pion masses) does not apply to pion 
production in few-nucleon systems, since the initial nucleons have a large relative momentum 



P = yjMm n + ^ ~ 360 MeV/c. Because of this an expansion in P/M = ^Jm n /M ~ 0.4 
was adopted Q|. But the convergence of this expansion is slow, and it is necessary to include 
several poorly constrained higher-order mechanisms in order to achieve a good description of 



the data, e.g., the short-range interactions and/or off-shell amplitudes mentioned above 21|. 
These effects are included in meson-exchange-model calculations of pp — > ppir , but even so 
such calculations struggle to describe spin observables 0, I22I . The present situation, with 
a detailed review of this interesting field, is given in Ref. |23j. As that review points out, 
in spite of its difficulties, the expansion in powers of ^Jm n /M does predict the correct 
hierarchy of mechanisms for the production of pions in a relative s-wave. It also successfully 
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explains p-wave pion production 

The np — > dir° reaction is of particular and current interest. The total cross section was 
measured some time ago at TRIUMF j^j], and fitted to the expression 

a{np^dir°) = ^(ar l + l3rj 3 ), (1) 

where rj = q/m^ is the reduced cm. pion momentum and Ref. 3| gives the values a = 
184 ± 5 /j,b and (3 = 781 ± 79 /xb for what they call the reduced (pion) s- and p-wave cross 
sections (but see also below). The factor 1/2 results from a comparison with the (Coulomb 
corrected) isospin partner pp — > dir + . The n + d — > pp data (applying detailed balance) are in 
agreement with np — > dir° [26]. However, more recent data on the pp — > dir + reaction gives 
a = 208±5/ib, [3 = 1220 ±100 fib Q and a = 205 ±9 fib, [3 = 791 ±79 fib Q, indicating 
a possible isospin violation. A model calculation (2^ ma nag ed to reproduce the data of 
Ref. 



f. Q, 

but it was subsequently pointed out that Ref. |29| overestimated contributions 
from short-range pion-production mechanisms [3(|. More recently the charge-symmetry- 
breaking (CSB) forward-backward asymmetry was measured to be = [17.2 ± 8(stat) ± 
5.5(sys)] x 10~ 4 (jy] . The asymmetry is an interference between charge independent (CI) 
and CSB s- and p-waves. Ignoring relative phases it can be expressed as 



Afb(np -> dn ) ~ + ^ ~ w _ v 4. Vj (2) 

an + prf V a a 

where a and (3 are the CSB equivalents of a and f3. It is believed that this asymmetry, 
together with the recent 5cr result of the dd — > an experiment j^, wm belp constrain 



some of the lesser known, but important, CSB parameters in the chiral Lagrangian |33l. l34j . 
For this reason there is a need for consistent chiral calculations for both reactions. A first 
report on the demanding dd — > air calculation has been published 34J, and a calculation 
incorporating some of the physics mandated by chiral symmetry that affects A^(np — > 
dii ) exists 33j, but a reliable extraction of CSB parameters necessitates a consistent chiral 
calculation of the reaction np — > dir . 

In this paper we will use the power counting appropriate for pion production to establish 
the leading- and next-to-leading-order (LO and NLO) diagrams that contribute to the CI 
s-wave np —> dir amplitude close to threshold. This is a preparatory, but important, step 
toward the calculation of the CSB asymmetry, since Eq. (J2J) makes it clear that without an 
understanding of the CI cross section the asymmetry is not under control. 

The organization of the paper is as follows. We introduce the theoretical framework in 
Sec. |n] The LO operator and results are presented in Sec. Hill while the NLO operator is 
introduced in Sec. IIVI The results at this order are discussed in Sec. EJ and the energy 
dependence of the s-wave cross section in Sec. IVII We conclude in Sec. IVHI 



II. THEORY 



A. General considerations 

The np — > dir reaction can easily be decomposed into partial waves. The lowest CI 
partial waves are 3 PiS, 1 SqP, and 1 D 2 p in the spectroscopic notation 2S+l Ljl n , where S, L, 
and J are the total spin, orbital, and total angular momentum of the incoming nucleon pair 
and Ztj. is the orbital angular momentum of the emerging pion. The lowest CSB partial waves 
are 1 P\s, 3 S\p, 3 Dip, and 3 D 2 p. 

As we will restrict ourselves to near-threshold CI s-wave pion production, only the 3 Pis 
transition will be evaluated in this paper. This wave contributes to a, while 1 SqP and 1 D 2 p 
both contribute to (3. Note, however, that a and (3 do not give an unambiguous separation 
of s- and p- waves. An s-wave amplitude that is dependent on the pion energy, e.g., the WT 
term, will produce a term quadratic in rj. This term's interference with the //-independent 
piece of the 3 Pi s-wave amplitude mimics the rf dependence of p-wave contributions to a. 
A separation of the rj 2 piece of the 3 Pi s-wave amplitude and the p-wave contributions can 
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be made using the energy dependence of the analyzing power A y . Defining the calculated 
s- and p-wave contributions a(n) and /3(n), their energy dependence can be written as 

a(n) = a + a'n 2 , 

Ptn) = Po + P'v 2 , (3) 

and the analyzing power as (again ignoring relative phases) 

A,, 



{a + {& + p )n 2 )n " \ 





M 


—v 


a 





2(3 2a 6lq J 



(4) 



This expansion is valid as long as j^n 2 <C 1, ~^n 2 <C 1, and a rf 1. (In fact, existing 
data for pp — > dir + are consistent with A y oc r) for r\ < 0.5 
With our normalization, the total np — > c?7r° cross section is 

?7 m n m p m n 1 



2 47r 2 s p 4 



E l-^l 2 > (5) 



spins 

where p is the relative cm. momentum of the incoming nucleons, -y/i is the total energy in 
the cm. system, and the sum is over the nucleon and deuteron spins. A general form for 
the np — > dn matrix element, Ai, in p-space, is then 

M = ^2E d J ^0f^ d (p'')A(p^ np ( P \ p), (6) 

where Ed is the deuteron energy, ipd{p") and ip n p(p',p) ar e the deuteron and np wave 
functions, p = p" — p', and A(p) is the pion-production amplitude. A recurring problem 
in any calculation of threshold pion-production amplitudes is the mismatch between initial- 
and final-state nucleon momenta. Since the initial loop momentum p' ~ p ~ P m 360 MeV, 
while the final loop momentum p" ~ 45 MeV, the two integrals in the matrix element need 
to be performed in quite different momentum domains. Critical questions are: Where is 
the large initial-state momentum absorbed in order to effect the transition between the two 
regimes? Does it have to pass through the deuteron wave function or are there pieces of the 
operator that can help absorb it? 

B. Chiral perturbation theory 

We will use the "hybrid" approach to derive the prediction for np —* dir° at thresh- 
old. In this approach chiral power counting is applied to the transition operator A. The 
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wave functions in Eq. (JBJ) account for the effects of initial- and final-state nucleons propa- 
gating nearly on their mass shell. Such states have small energy denominators that do not 
obey the usual ;\RT power counting j^]. Thus they have to be calculated non-perturbatively, 
ideally in a manner consistent with chiral symmetry. The hybrid approach has been applied 
successfull y to a wide range of pionic and electroweak reactions in the two-nucleon system 



'see Refs. 136 
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for reviews). 

We use the chiral Lagrangian as developed in Ref. |17| . Eqs. (3) and (6). These contain 
pieces from the first and second chiral orders, which is all we need for the present purposes. 
The chiral Lagrangian is used to construct the transition operator A in Eq. (jHJ, which is 
just the sum of all of the two-nucleon-irreducible diagrams that contribute to NN — > NNir. 
These diagrams are arranged using a power counting which involves a number of different 
energy /momentum scales. Ordered by size these scales are: 

• The pion momentum |q| which is <C m n close to threshold. At NLO (which is all we 
consider here), this scale does not enter the s-wave production amplitude. 



The deuteron binding momentum 7 = \J BM = 45 MeV, where B is the deuteron 
binding energy. We expect that the integral over p" in Eq. @ will be dominated by 
|p"| ~ 7. However, here we will conservatively adopt |p"| ~ m n . 

The pion mass m n for which we take 138.039 MeV. 

The Delta-isobar-nucleon mass difference A = Ma — M = 290 MeV ~ 2m 7r . 

The nucleon momentum needed to produce a pion at threshold P = 



yJmjM + m^/A) = 360 MeV. 

The chiral symmetry breaking scale A x ~ 1 GeV. 



We will see below that the scale A will never appear explicitly in the tree-level diagrams, 
since it will be combined with the pion mass in propagators. Meanwhile, in loop diagrams 
involving the Delta-isobar, we follow Ref. 3^| and assign A ~ P. 

Thus we have a three-scale problem, with the hierarchy m n C P < A x . We choose to 
expand our amplitudes in the dimensionless parameter 

^ = = \ tt ~ ~; — ~ 0.4. (7) 
P V M A. w 



When counting powers of x in two-nucleon-irreducible loop diagrams it is important to 
distinguish between graphs of the type in Fig. and those of the type shown in Fig. El In 
the first case the large nucleon momentum P may pass through the nucleons and exchanged 
pions only (Fig. [TJ, whereas in the second it must go through pions in a loop [e.g., Fig. Efa)- 
(h)]. In consequence the loops in Fig. [I] are independent of the initial nucleon momentum P, 



FIG. 1: Examples of diagrams where the large momentum need not appear in a pion loop. 

since the loop resides only on a nucleon which is moving with momentum P. By translation 
invariance, the loop integral then scales with m n only. In Fig. |3] the loop integrals pick up 
the large momentum scale and the integral will be dominated by pion modes of momentum 
P. 



C. Wave functions consistent with chiral symmetry 

In evaluating the matrix element © we would like to employ wave functions which are 
consistent with the chiral power countin g th at we use for the operators. One way to do this 
for the deuteron was developed in Ref. 40(. This method was extended to NN scattering 
in 4l| . In these wave functions the long-range part of the NN potential is derived using 
saa. At the same time our ignorance about the short-distance physics of the 
nucleon-nucleon interaction is parametrized by matching, at some chosen radius R, the 
wave function obtained for r > R to a square well solution for r < R. Here we will use this 
approach to construct the initial-state ny wave function, starting from the 3 Pi np phase shift 



of the Nijmegen phase shift analysis [42j. Since we are dealing with cm. energies of order 
m w , we determine the asymptotic scattering wave functions from the Nijmegen phase shifts 
at a lab. energy of 279.5 MeV and then expand around this point. This is different from the 



procedure chosen in |41| . where, due to the low NN energy, the nn scattering length and 
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effective range could be used as input. Note that, at the present state of the calculation, 



we include only one-pion exchange (OPE) with f% 



NN 



0.075 



43| as the long-range part 



NN potential. In the future we will also implement chiral two-pion exchanges 
. (Unfortunately, it is impossible to use high-order chiral potentials directly, 
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47| to 



of the chiral 
(TPE) 

since in general they are accurate only for energies below 250 MeV |7|,|45j.) As a cross check, 
below we also use some of the available low-% 2 potential model wave functions 
calculate the matrix element (JHJ). 

The chiral wave functions we have constructed for use in our calculation have a co- 
ordinate space regulator present in them: the square well of radius R. Other regulators 
(e.g., Gaussians or sharp cutoffs in momentum space) have been used in obtaining NN 



0,0,0,0,0- 



In all of 



wave functions from a potential derived via a expansion 
these studies a regulator is introduced before the potential is iterated using the Schrodinger 
equation. This means that when, in Eq. ©, matrix elements of pion-production operators 
are computed using such wave functions, there is a scale-dependent regularization in the p' 
and p" integrations. The scale dependence of this regularization is less obvious when NN 
potential models are used to calculate the t/Vs appearing in Eq. (JHJ). However, even there 
it seems clear that the hadronic form factors in such models somehow implicitly provide a 
regularization that is explicitly present when wave functions derived from the effective field 
theory are used. 



III. LEADING-ORDER OPERATOR AND RESULT 



The LO [or 0(x)} s-wave operator is given by the Weinberg- Tomozawa (WT) and the 
(nucleon recoil) one-body term; 



•Awt(p) 
Ab(p) 



3g A uj 



e r t T 2 



(T 1 ■ (p 



tr 2 ■ (p + 



-ig A uj q (2ir v ' 
4Mf w 



m; 2 + (p-a)2 



m; 2 + (p + a)2 



(p + |) a x ■ (p" + P ' - I 



fp-f)<r 2 -(p" + p' + f 



(8) 



where m'J 2 = \{rn^ — q 2 ). [In order to derive Eq. (JSJ) we have assumed equal energy sharing 
between the two nucleons.] 

Chiral power counting for the operator A will be a useful guide to the size of physical 
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mechanisms only if different contributions to the operator are sensitive to similar momenta in 
the NN wave functions used in the evaluation of the matrix element Q. This is clearly not 
the case for the operators in Eq. (jSJ): the one-body operator "lives" off the high- momentum 
components in the incoming state, while the two-body nature of the pion-rescattering mech- 
anism allows better momentum sharing, and access to the larger low-momentum components 
in the np and d wave functions. Therefore in order to estimate the impact the operator Am 



17J] and use the Schrodinger equation for the deuteron 



has on observables we follow Ref. 
state \ipd) f° write 

(^IAbIW) = (MVGoA 1B \lpnp)- (9) 

(Here V is the NN potential used to generate the wave functions.) Once Eq. © is invoked 
momentum sharing can take place through V. So, in order to power count the operator 
Aib, the pion exchange and contact interaction that make up V at leading order in x?T 
are included in the diagrams. This also puts the nucleons in the "external" states closer to 
being on-shell, and facilitates attaching wave functions that are calculated assuming on-shell 
nucleons. 



(a) 



(b) 



(c) 



FIG. 2: Leading-order diagrams for s-wave pion production on two nucleons. 



The diagrams in Fig.|2]are then the same order in chiral power counting, which justifies the 
statement made at the beginning of this section that Aib and -Awt together give the leading- 
order amplitude. This statement is borne out if we compute the different contributions to 
«o from these leadin g-or der processes using a variety of phenomenological and chiral wave 



functions 



40 



46 



47 



5ll | . The results of that computation are shown in Table HJ 
Considering that the power counting predicts the WT and one-body term to be of the 
same order, the six order of magnitude difference between the two in the first row of Table H] 
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TABLE I: Values for ao obtained from various different contributions to the leading-order np — > 
amplitude. All results are in units of fib. Note that in each case (except the last column) the result 
for ceo is the one obtained by considering just that mechanism. 



potential 


ao(WT) 


ao(lB, PWIA) a (lB, ISI) 


a (lB) 


«o(LO) 


Nijm I 


124.0 


11.94 


11.83 


0.0003 


123.6 


Nijm II 


96.8 


16.26 


18.21 


0.0682 


102.0 


avl8 


82.2 


15.75 


22.35 


0.6961 


98.1 


Bonn A 


111.8 


6.264 


9.622 


0.4234 


126.0 


Bonn B 


110.4 


11.28 


12.96 


0.0713 


116.0 


LO xPT (R = 1.5 fm) 


40.6 


5.809 


19.90 


4.723 


73.0 


LO xPT {R = 2 fm) 


114.7 


1.124 


9.391 


4.325 


163.6 



should cause some concern. However, we can split up the matrix element of Am into a 
plane-wave (PWIA) and initial-state-interaction (ISI) part, i.e., write: 

(ipd\A 1B \ipnp) = (MAib\p) + (MAmGoVtyr*). (10) 

The fact that the PWIA and ISI contributions to do are of the same size lends credence to 
our application of power counting to the pion-production operator. The one-body ISI and 
PWIA values alone are closer to the do from the WT mechanism than is their sum, although 
each piece still differs from a (WT) by an order of magnitude. The extremely small matrix 
element on the left-hand side of Eq. (|10|) apparently results from a cancellation between two 
numbers of a little smaller than natural size on the right-hand side of that equation. The 
total LO value for do is then quite close to the number from Fig. Efa) alone. 

The numbers in Table H] are in accord with those given in the earlier study of da Rocha 
et al. 0| . (Cancellations between different pieces of the integral, leading to a reduction 
in the one-body term, were also observed in Ref. Q , although there the PWIA and ISI 
contributions are not given separately.) 

We find that there is considerable wave-function dependence in the result for olq. The 
calculation with the chiral wave function with cutoff R = 1.5 fm is expected to differ from 
the other results, since at p — 360 MeV the matrix element will probe distances where the 
neglected TPEs become important. If we employ only the low-^ 2 NN potential models we 
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find a moderate spread in values: 

a = 111 ± 13 fib. (11) 

This variation with the choice of NN model reflects the sensitivity to short-distance physics 
in the LO pion-production matrix element, since these potentials reproduce the 1993 NN 
database with \ 2 ~ 1 and have the same OPE tail. The LO prediction (jllj) gives « 60% 
of the measured a, which is consistent with the accuracy expected from a LO result in an 
expansion in ^Jm n /M. 

The use of wave functions derived assuming zero-energy transfer between the nucleons 
means that there is an inconsistency between these wave functions and the operator. In 
diagram |2^b), for example, we assume there is energy transfer between the nucleons, as 



dictated by the Feynman rules and our assumption o 



equal-energy sharing between the 



3 



for power counting diagrams by 



initial particles. If the procedure used here and in Ref. 
"pulling" the last pion exchange from the wave functions is accurate then this inconsistency 
should be a higher-order effect, since it involves the difference 

1 1 = m l ~ < 2 x 1 (12 ) 

m; 2 + (p±|) 2 m^ + Cpif) 2 m^ 2 + (p±§) 2 m ; 2 + (p±|)2- ^ > 

The momentum transferred here is of order P, so the difference between including and not 

including energy transfer in the one-pion-exchange propagator is of relative order , as 

long as q ~ m v . It therefore contributes first at two orders beyond leading, which is one 

order higher than we consider here. 

Similar arguments ensure that any rj dependence of the leading-order s-wave production 

amplitude A s arising from energy flowing through these meson denominators is also a NNLO 

effect. Thus 

a, = uj q [A, + e>0V)] > (is) 

where A s is the reduced leading order As and the 0(x 2 r] 2 ) correction stems from the q 
dependence of the vertices and propagators in Eq. (JBJ). Thus the prediction of %PT for these 
diagrams ()13jl is that (ignoring ISI for the moment) 

a = ^a [l + r ] 2 + O(x 2 r ] % (14) 

which implies that a' = 6lq at LO. However, a numerical calculation including the full 
scattering state reveals that 

a' = -57±4/ib (15) 
11 



for the potential models, and somewhat smaller, but positive (a! = 41 ±30 fib), for the chiral 
wave functions. The energy dependence coming from ISI apparently provides a > 100% 
correction to Eq. (|14j) . We will return to this issue of energy dependence once the NLO 
contributions have been calculated. 



IV. NEXT-TO-LEADING ORDER OPERATOR 

At NLO [or 0(x 2 )}, there are pion loops and a Delta-isobar-excitation term. The cor- 
responding diagrams are shown in Fig. El Some of these were discussed in Ref. j^2|, and 
a complete treatment can be found in Ref. j^J. Here we draw on the analysis of Hanhart 
and Kaiser for the pion loops shown in Fig. since for these graphs they employed the 
power counting in x we are using here. 

At NLO and in threshold kinematics diagrams |3fb) andEfc) cancel for NN — ■> dir, while 
Efe) andlH^f) oc m^-lnm,,-, and thus are of higher order [3,9] • Diagrams Efg) and^h) are 
also of higher order in the x expansion because of their vertex structure. As for diagram 
Efi), since the deuteron is isoscalar this particular Delta-isobar diagram cannot contribute 
to NN -> dir. 

There are also a number of Delta-isobar loop diagrams which nominally appear at NLO, 
but they cancel each other in threshold kinematics. This is associated with the absence 
of a chirally-invariant counterterm at this order j^J. Thus for s-wave pion production in 
np —> dir°, there are only two loop graphs that need to be calculated if an answer accurate 
to 0(x 2 ) is desired. 

The sum of these two diagrams (a) and (d) can be evaluated straightforwardly. The 
corresponding amplitude, together with that for the analogous diagrams in which the roles 
of nucleons 1 and 2 are interchanged, gives the following result for the operator mediating 
the NN — > dn transition: 

Aa)+(d)(p) = \{<Tl + ^ • P« - O^HpT), (16) 

with gA the axial coupling of the nucleon, f n the pion decay constant, and B given by the 
scalar integral 

•~2n _ 9a 1 f d4 Q (q + P) -q 



B{P ) " 8/6 i J (2vr) 4 vq[(p + qf - ml] [q 2 - ml] ' 
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(a) 



(b) 



(c) 



(d) 
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(h) 



(i) 



FIG. 3: Next-to-leading-order diagrams (with no Delta-isobars in loops) for s-wave pion production 
on two nucleons. 

Here v is the four-velocity of the heavy fermions in the non-relativistic %PT Lagrangian. 
Since v ■ p ~ m n for threshold pion production, to the order to which we work B can be 
regarded function of p 2 alone. 

If that function is computed using dimensional regularization (DR) we find: 



BMP) 



-p 2 + 0(m n ) 



(18) 



Clearly if p 2 = —Mm^, as is the case in threshold kinematics, this expression vanishes in 
the limit — > 0. This is consistent with the chiral-symmetry requirement that the pion 
should couple softly to the iVTV system. 
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However, to get the result (|18|) the fact that the scale-less and linearly divergent integral 

d d l 1 



(19) 



(27r) d (v ■ l)P 

is zero in DR has been employed. If a scale-dependent regularization scheme is used this 
integral will no longer vanish. For instance, with Pauli-Villars regularization we find 



256/s 



4A 

7T 



-p 2 



(20) 



which does not vanish in threshold kinematics as the chiral limit is taken. We will return to 
this point at the end of Section IV! 

The result (II 6 j) and ()18j) agrees with that given in Ref. |29j if exchange diagrams are 
added to the computation. Here, instead of calculating such diagrams explicitly, we include 
the effects of fermion anti-symmetry by only summing over those initial- and final-state iViV 



partial waves which are allowed by the Pauli exclusion principle (see Section III Ajl . The 
result (llfjjl for the operator A is thus the correct one to employ in Eq. (|5jl. since our A^A" 
wave functions are calculated on the (normalized) isospin basis. 



V. NEXT-TO-LEADING ORDER CALCULATION 

After Fourier transforming to r-space the matrix element flOJ) becomes: 

M = ^2E d J £r^ d (v)A(r)iP np {v, p), (21) 

where A(r) = J -^^e l ^' r A(p). The calculation of this Fourier transform for the s-wave 
tree-level diagrams is straightforward, but care is required when the pion loops are included. 
The Fourier transform of Eq. (JTfiJl can only be carried out after a p-space regularization [e.g., 
by including a factor exp(— p 2 /2A 2 ) and then letting A — > oo]. The result is 

Aa) + (d)(r) = \{<r l + a 2 ) ■ r« " ^Sf^ (22) 

This highly singular operator causes a linear divergence in the matrix element ()21j) . as can 
be seen from the following arguments, which are valid for all wave functions calculated from 
a local AW potential. At small enough r (< i?div, say) the radial deuteron s-state wave 
function u(r) ~ kq^, while the radial 3 Pi np wave function v\(r) m K\r 2 '. (Here the Ki 
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are constants and i?div indicates where these approximate forms no longer hold.) Thus the 
short-range contribution to the matrix element behaves as 

f Rdiv dr 

kqKi / — , (23) 
jo r z 

which is clearly divergent. Thus the matrix element [in contrast to the operator Eq. (|16j)] 
does not vanish in the chiral limit, since the wave functions will still have the above behavior 
as r — > even if the pion-production threshold becomes E = (though we would have a 
different value for K\ in that case). The loop matrix element violates chiral symmetry! 

In order to render the integral (|2*3*|) finite we introduce a new (coordinate-space) cutoff 
regulator 1 — exp(— ( 2 r 2 ). We find 

r dr <r)v{r)[l-eM-e^] ^ + Q 

Jo r 5 

for large enough (, since the integral is then dominated by its behavior at small r. The 
constant C is wave-function dependent, but independent of ( as long as -RdivC ^ 1- The 
linear cutoff dependence predicted in Eq. (|2*^|l is confirmed in the actual calculations using 
a variety of wave functions and cutoffs (, see Fig. @] 

In the context of pionic and electromagnetic processes in the two-nucleon system the 
pion-production reactions we are considering here are special. In all of the other reactions 
that we are aware of the operators occurring at low order in the chiral expansion diverge 
much less rapidly as r — > than the \ of Eq. (J23). In consequence the matrix elements 
of these operators do not have the sensitivity to the short-distance piece of the NN wave 
function seen in Fig. 01 For an elegant demonstration of this in the context of pion-deuteron 
scattering, see Ref. [531 ] . 

The linear divergence of Eq. (|2*^j) (or Fig. is removed by introducing a short- distance 
operator: 

1 / I \ ~t-a • a\^o(0 
+ tr 2 ) ■ pfa - %T 2)-^i 

where g?o(C) i s a dimensionless, cutoff-dependent constant. The operator (|23|l does not come 
from any term present in the original Lagrangian of but such a contribution is required 
here because cutoff regularization is employed to compute the matrix element (|21j) . 



At(p) ^^1 + <T 2 ) • P(t7f - iT%)-^, (25) 



j3 



The operator (|25|) was not included in the analysis of Cohen et al. 
does not couple derivatively (or via powers of the pion mass) in A c t, and there 
operator cannot be derived from the Lagrangian of The analysis of Ref. 

15 



because the pion 
ore this 



17| uses 



10 



8 - 



•K 6 

a 

s 

hi 
.a 
- 



g 4- 



2 - 



1 1 1 1 1 

— R = 2.0fm 


■ 


1 1 1 


1/ 


- R = 1.5 lm 








— — Nijm I 






— 


-- Nijm II 








.... qvI a 






y 






y 

y 

y 

y 

y 


y 

y 

y 

y — 

y 


- y 


' / 
y 

y 

y 

y 

y 




— 
- 


/ y 
/ y 








/ y 
— X y 
/ y 








/ y 
/ y 
/ y 
/ y _ • — 

/ y — • ~~ ' 


















. 







5000 10000 
cutoff (Me V) 



15000 



FIG. 4: Cutoff (() dependence of the loop matrix element (|21j) for various wave functions as 
indicated. Two chiral wave functions with different matching radii R were used, as well as some 
low-x 2 potential models. 

that Lagrangian, together with naive dimensional analysis (NDA), to order all of the pion- 
production operators that are consistent with QCD's chiral symmetry and the pattern of 
its breaking. In this case the regulator we have used when Fourier transforming the loop 
am plit ude (J16|) breaks the chiral symmetry invoked when NDA was applied by Cohen et 
al. |l2(. Therefore we expect to have to add extra operators which break the symmetry in 
the same way as the regulator. Including such terms allows us to write down a renormal- 
ized amplitude that respects the symmetry — even though neither the loops or counterterm 
individually do. The term f)25j) must be included in our calculation, in spite of the fact that 
it violates chiral symmetry, so that it can cancel the chiral-symmetry-violating divergence 
(J2IJ) that is proportional to the cutoff (. 

Calculating the matrix element of A c t in r-space, we find [via Eq. (12 lj) ] 

6d (C) 



—el 



IP x e npJ 
1(3 



2E d K K 1 , 



(26) 



where e np is the polarization vector of the spin-1 nucleon pair and the Fourier transform 
creates a <5^ 3 ^(r) which picks up properties of the wave functions at the origin only. This can 
be evaluated without any regularization. It should come as no surprise that the expression 
(J26j) also does not vanish in the chiral limit. Of course, one can always choose the value of 
do such that the sum of matrix elements for the loops and the counterterm vanishes in the 
limit m n — > 0. However, since we are interested in pions with observed masses, we choose 
to extract c?o(C) & t the physical threshold. This is to say, we fix the value of c?o(C) by fitting 
the experimental value of the threshold cross section |25j. Given the behavior of the loop 
integral, do will then run linearly with the cutoff (. By construction the total matrix element 



will now reproduce the value of a from Ref. 



. Thus the renormalization not only removes 



the divergence in the np — > dir° matrix element that occurs at NLO, it also removes the 
wave-function dependence that we found in the LO result in Section IIHI 

An alternative regularization of the matrix element of graphs a) and Efd) is provided 
by refraining from taking A — > oo directly after Fourier transforming, and instead keeping 
the regulator mass finite until after the matrix-element integral (f2"T]) has been performed. If 
we do this with a Gaussian regulator exp(— p 2 /2A 2 ) the loop matrix element becomes 



-Mioops = -4 • ( P x e np )^- 5 ^A 5 J dru(r)I(Ar) Vl (r), (27) 

where I(x) = ±[± - x + y/2(£ + 2- x 2 )F(-x/V2)) and F(x)[= -F(-x)] is the Dawson 
integral 



2 / X . j-2 



F(x) = e~ x / dte 1 . (28) 
Jo 

Surprisingly, in this procedure the coefficient of the linear divergence vanishes for local wave 
functions. But, it is still the case that the value found for the loop matrix element is sensitive 
to the short-distance physics. The result is strongly dependent on the wave function (see 

Fig. ED. 

The term linear in A vanishes because there is a cancellation between the short- and 
long-distance (as separated by 1/A) behavior of the integrand. If a function F^(p) is used 
to regulate the Fourier transform then the divergent piece of the integral is [for fixed A and 
using the r < R div limits of u{r) and t>i(r)]: 

Raw d f d 3 p - 

dru(r)vi(r) — / ———pe ip ' r F^(p) 



drJ (2tt) 3J 
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FIG. 5: Single-regulator (Gaussian) cutoff (A) dependence of the loop matrix element for various 
wave functions as indicated. 



where we have taken the limit A — > oo in the second term on the last line, and F\(r) = 
J °° dpsin(pr)F\{p). Any regulator for which the square bracket in the last line vanishes (a 
class which includes exponential and Gaussian regulators) produces a A-independent result 
in the A — > oo limit. Therefore the coefficient of the linear divergence is identically zero 
for such regulators. And indeed, Fig. |5] (calculated with a Gaussian cutoff) shows that 
calculations with a variety of wave functions produce flat behavior for 1-2 GeV ^$ A < 
10 GeV. (Qualitatively similar curves result if an exponential cutoff is used.) 

For the local potentials, the flatness region persists at least up to ~ 200 GeV. Non-local 
potentials (e.g., Nijm I, Bonn models) start to deviate from flatness around 10 GeV, a hint of 




J r=R di 



(29) 
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which is visible in Fig. |S] Their non-locality introduces additional high-momentum compo- 
nents that influences the wave functions at small r. These high-momentum components are 
neglected in the two-step procedure and there is no distinction between local and non-local 
potentials in Fig. EJ 

The results of Fig. El confirm the prediction of Eq. (|29|). They do not show a divergence as 
the cutoff mass A is taken to infinity. However, the matrix element remains non-vanishing in 
the chiral limit, and its value for m n = will be determined by the high-energy mass scales 
k , K\, and -Rdiv These mass scales are related to the (implicit or explicit) scale-dependent 
regularizations carried out in obtaining the NN wave functions we have used to calculate 
the matrix element. In order to remove this dependence on the (unphysical) wave-function 
regularization mass scale(s) we again introduce a short-distance operator of the structure 
f!25|) . But this time we retain the regulator after Fourier transforming the operator to co- 
ordinate space, treating it in the same way as we did the loop-diagram operator when we 
obtained f)27j) . The matrix element of the short- distance operator then reads: 

M ct = -2eJ ■ (p x e np )^-y^A 5 [ rdru{r)e-^ 2 '\ h {r). (30) 

We now include the short- distance operator (J3U|) and adjust do(A) such that the total 
matrix element squared gives «o = 184 pb. This yields the values for do(A) shown in Table ITT1 
The (io's shown there are, for most wave functions, quite rapidly varying with A. This is 
because the values given are for A below the point where -Mi 00 p becomes A independent, see 
Fig.0 

For each A there are two solutions for do since we are solving a quadratic equation. Once 
p-waves are included in the calculation it should be possible to resolve this ambiguity with 
polarized data, e.g., A y , which is available for pp — > dn + at low energies j^]. Such data 
fixes the relative phase between the s- and p-waves, provided the phase of the p-wave is 
known, see Eq. (jlj). We reiterate that, regardless of which value of do turns out to be 
correct, the renormalization that is necessary at NLO in the threshold np — > dn amplitude 
automatically eliminates the wave-function dependence that is present in the LO amplitude. 

This procedure of regularizing using a Gaussian and keeping A finite until after the 
matrix element is calculated seems more intuitively appealing and consistent than the "two- 
step" procedure introduced above. The fact that the term linearly proportional to the cutoff 
vanishes in such a "one-step" regularization also means that renormalization does not involve 
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TABLE II: Values for do for various wave functions and cutoffs, using "one-step" regularization 
with Gaussian. For each A there are two solutions, labeled (1) and (2). 





A = 400 MeV 


A = 800 MeV 


A = 1200 MeV 


JJVJ IjCll Lldl 


A (A \ 




«o(iJ 


A / r >\ 

do{z) 


A (-\\ 

«o(,ij 


A 1^ 


l\Jl 1 TTI I 


1 R9/I 


n DOR*:! 


1.00 I 


n 19c; 
-u. 1ZO 


1 .Out* 


n 1 fin 

-U. 1DU 


Nijm II 


1.742 


-0.199 


2.629 


-0.482 


4.780 


-1.079 


avl8 


1.928 


-0.2536 


4.139 


-0.857 


10.138 


-2.320 


Bonn A 


1.697 


-0.1017 


1.944 


-0.194 


2.365 


-0.291 


Bonn B 


1.670 


-0.1325 


2.006 


-0.257 


2.678 


-0.437 


LO xPT(1.5 fm) 


3.331 


-0.7221 


4.103 


-0.935 


4.368 


-1.004 


LO x?T(2 fm) 


1.887 


0.0072 


1.502 


0.0169 


1.445 


0.0179 



taking the difference of two large numbers. 

Nevertheless, there remains a concern regarding the consistency of our use of two differ- 
ent regularization schemes: DR to compute the operators and cutoffs for the integrals we 
evaluate to obtain the matrix elements of those operators. Ideally the same regularization 
should be used for both pieces of our calculation. In this context we remind the reader that 
introducing a Pauli-Villars regulator (instead of DR) in the evaluation of the diagrams Eta,) 
and (d) results in a linear divergence appearing in the function B(p 2 ) [see Eq. ()20j) ]. There- 
fore if a scale-dependent regulator is used to compute the loops of Fig. El a counterterm of 
the form (|25|) must be introduced already at the level of the operator, in order to facilitate 
renormalization of this linear divergence. Using DR — in which the renormalization of such 
linear divergences is performed implicitly by setting scale-less, linearly divergent integrals to 
zero — to compute B(p 2 ) postpones the necessity to introduce this counterterm. But it can- 
not be prolonged indefinitely if scale-dependent regulators are used to calculate the matrix 
element of the operator (|lfijl . 

VI. ENERGY DEPENDENCE OF a AT NLO 

The NLO loops are proportional to — p 2 ~ Muo q . Therefore, as at LO, xPT predicts 
that a! = a>o for the PWIA piece of the amplitude. The advantage at NLO is that short- 
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range pion-production operators are included, which helps ameliorate the mismatch between 
the initial- and final-state momenta. Of course, a is now correctly reproduced, since the 
counterterm values given in Table ITD are found by adjusting a to the experimental 184 /xb. 
Calculating the energy dependence of the reduced s-wave cross section and fitting it to 
a = 6lq + a/77 2 for various A's, gives the values of Table ITTT1 (For cutoffs above ~ 1 GeV the 
calculation becomes independent of A.) 

TABLE III: Values for a' (in /A>) for various potentials, using counterterm values extracted for 
different A. The calculation is done for both do solutions (1) and (2). 





A = 


400 MeV 


A = * 


100 MeV 


A = 


1200 MeV 


potential 


a'(l] 


a'(2) 


6/(1) 


a'(2) 


a'(l) 


a'(2) 


Nijm I 


29 


-61 


98 


-53 


108 


-52 


Nijm II 


54 


-78 


144 


-62 


169 


-59 


avl8 


28 


-84 


102 


-72 


121 


-69 


Bonn A 


14 


-66 


78 


-59 


87 


-57 


Bonn B 


20 


-72 


82 


-64 


92 


-63 


R = 1.5 fm 


777 


205 


1471 


352 


1637 


383 


R = 2.0 fm 


277 


93 


434 


103 


466 


104 



Regardless of which choice we make for do, the chiral wave functions give results for a' 
that deviate considerably from the potential models and from each other, which might seem 
to be a disconcerting result. However, this deviation presumably comes from our neglecting 
TPEs, which seem to play an important role in processes at this momentum transfer. Our 
chiral wave functions, especially the one with R = 1.5 fm, give a poor description of physics 
at r < 2 fm. With this in mind we conclude that the range of NLO predictions for the two 
solutions is: 

a (1) = 149 ±41 /ib, (31) 
6/(2) = -60 ± 8 /ib. (32) 

The fact that (in either case) a' 7^ oo can be attributed almost entirely to distortions in the 
initial state, since higher-order q dependence in the amplitudes is a much smaller effect. 

In fact, using the first d value yields a's in reasonable accord with the PWIA ansatz 
a' ~ a . If, on the other hand, the second solution is chosen, the NLO prediction for a' 
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has only a very small shift from the LO result (|15j). Because of this convergence, we favor 
the second solution. The chiral series for a! seems to be behaving well, so we anticipate 
that higher-order corrections to the prediction (|32|) are < 15%. Once the p-waves have been 
calculated this result can be tested by comparison with, e.g., the analyzing power. 

In any case, the results of Table 11111 indicate that some of the experimental /3 does arise 
from the energy dependence of the s-wave pion-production amplitude. However, p-wave 
pion production gives the greater part of the observed number. 

VII. CONCLUSIONS 

We have calculated the s-wave cross section for np — > dir (or, equivalently, the Coulomb- 
corrected pp — > dir + ), using chiral perturbation theory to NLO. We find that the dimension- 
ally regularized NLO pion loops give a divergent result when they are sandwiched between 
wave functions, violating chiral-symmetry constraints. Thus the matrix element itself needs 
to be regularized. The variation with cutoff of this matrix element depends on the choice 
of regularization procedure. Regardless of the particular procedure chosen, a short-distance 
operator is needed to renormalize the theory. Matrix elements of both the loops alone and 
the short-distance operator alone are dependent on short-range physics, i.e., they display a 
strong variation with choice of wave function. However, by construction, their sum is cutoff 
independent at the pion-production threshold. A similar sensitivity to short-range physics 
has been observed in a calculation of pp —> ppn j^j]. 

Therefore it is not immediately apparent that ;\RT can make more than a leading-order 
prediction for the threshold amplitude in np — > dir . This prediction has a modest wave- 
function dependence, but, being a LO result, is only accurate at the ±40% level. However, 
%PT does predict the energy-dependence of the s-wave pion-production operator. With 
an appropriate choice of the counterterm, this prediction is very little changed by NLO 
corrections. Pions produced in a relative p-wave are expected to give the remainder of 
the (non-trivial) energy-dependence of the near-threshold np — > dir cross section. Under- 
standing charge-independent p-wave pion production near threshold is also critical to an 
interpretation of the TRIUMF result for A^(np — > dir ), and so we are currently working 
on a calculation of the amplitude for p-wave pion production in np — > drr . 

Note: Immediately before submitting this paper we became aware of an ongoing study by 
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Lensky et al. 3]. They propose a solution to the difficulty with the NLO pion-production 
operator we have identified here. This solution involves a careful treatment of the two- 
nucleon-reducible and two-nucleon-irreducible pieces obtained when the leading-order mech- 
anism depicted in Fig. |2[a) is sandwiched between NN wave functions. Presumably this 
also leads to a different prediction for the energy dependence of the NLO s-wave np — > dn° 
cross section. 
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